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We study a dynamical system model of an idiotypic immune network with a small number 
of degrees of freedom, mainly focusing on the effect of a threshold above which antibodies 
can recognise antibodies. The response of the system to invasions by antigens is investigated 
in the both models with and without the threshold and it turns out that the system changes 
in a desirable direction for moderate magnitude of perturbation. Also, the propagation of 
disturbance by an antigen is investigated in the system of one-dimensionally connected basic 
units taking the closed 3-clone system as a unit, and it is clarified that the threshold of 
the system has effects to enhance the stability of the network and to localise the immune 
response. 

Introduction 

It is experimentally well-known that an immune system is activated by interact- 
ing with itself to prepare for unknown antigensEP'tf. Takingthis fact into account, 
an immune network model was introduced by N. K. JerneEP. Later, a dynamical 
system model of the immune network was introduced by F. J. Varela et. al.&. In 
this model, the basic elements are antibodies and B-cells which produce antibodies. 
The effects of helper and suppressor T-cells are taken into account by introducing 
functions which represent interactions between antibodies, and between B-cells and 
antibodies. This has been called the "second generation immune network model" 
In this model, some imuoriant_mechanisms such as recognition, memory, and toler- 
^ ■ ance have been studied tP < Er LP <u . 

In the previous paper B*, we studied the original Varela model and also the 
modified model. Mainly, we have investigated the model with a threshold above 
which antibodies can recognise antibodies and have reported that the system has 
chaotic states, and also a peculiar type of limit cycle, which we called differentiating 
state. Further, as for the response to the invasion by antigens, we have found that 
when the system is in a differentiating state, its response to a specific type of antigen 
is sensitive. That is, differentiating states are considered to represent a kind of short 
term memory of the invasion by antigens. Further, we have found that the response 
time in a chaotic state takes an intermediate value compared to differentiating states. 

An important issue which has not yet been addressed is the question of how 
central the use of thresholds is to the results found. This is investigated in the first 
part of the present paper. In the second part of the paper we go on to study the 
localisation of immune responses. It has been pointed out by De Boer et.al. that a 
fault of the network view of immune systems is that when only one clone interacts 
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with one antigen, immune responses spread over the whole system and the activation 
by the invasion by antigens is not localised in the model by Varela et.al. This be- 
haviour seems inappropriate because it suggests that large networks, such as might 
be found in the real immune system, can be unstable since small local perturbations 
can generate global effects by spreading over the whole system. Thus, it is very 
interestingto find a mechanism to remove this fault and to localise the immune 
responses©. Thresholds may be one of such candidate. Therefore, we study the 
effect of thresholds on the response to the invasion by antigens, by investigating the 
propagation of disturbances in one-dimensionally-connected basic units. 



Model 

The evolution in time for the concentrations of an antibody, /j, and of a B-cell, 
bi, with the idiotype i(i = 1, 2, • • • , N) are given by the following equations. 



dfi/dt 
dbi I dt 



-Kxuih - K 2 fi + K 3 Mat ( fff ) b, t 
-K 4 h + K 5 Prol 0*) bi + K 6 . 



(0-1) 
(0-2) 




See Ref. 9) for detailed definitions and 
explanations. Here, is the sensitiv- 
ity of the network to the ith idiotype, 
defined as 



N 



(0-3) 



Fig. 1. Profiles of Mat(a) and Prol(a) adopted 
in this paper 



The interaction between two different 
idiotypes i and j are represented by the 
connectivity m°-. We set > if 
there is an affinity between i and j, and 
m?- = if there is not. For simplicity, 



we assume m 



o 



jj, where <5$j is the Kronecker's delta. Mat (o~i) and Prol (o~i 



are the interaction functions between antibodies, and antibodies and B-cells, respec- 
tively. In general, Mat (crj) and Prol (o"j) are assumed to be unimodal functions and 
Prol (dj) to be shifted to the right with respect to Mat (cij) in order to take the roles 
of T-cells into account. The functions we adopt here are the followings. See Fig.l. 



Mat (o-j) 
Prol (<7j) 



Ui [tanh {U 2 {oi 
U4 [tanh {U5 (uj 



T lm )} - tanh {U 3 (o-i - T um )}] , (0-4) 
T lp )} - tanh {U 6 (o-i - T up )}} . (0-5) 



Model with Threshold 

First, let us study a 3-clone closed chain system with a threshold kq above which 
the ith antibody can recognise other antibodies. Each element of the connectivity 
matrix M = {rriij} is defined as 



mij(t) = m° ij O(f j {t) - kq), 



(0-6) 



where &(x) is the Heaviside function (i.e.,6>(x 
In this system, there exist differentiating states 



; 1 for x > and for x < 0). 
We consider the response of the 



Letters 



3 



system in differentiating states to invasion by antigens. In these states there exist a 
clone with a long period and a clone with a short period. We denote the former by 
L, and the latter by S, respectively. We treat two cases of invasion by antigens in 
the followings. 
Case 1 

An antigen A similar to the antibody fa invades the system. We assume in this 
case that antibodies fa and fa recognise the antigen as the antibody fa. Thus, in 
the sensitivities 02 and 03, f\ is replaced by f[ = fa + A. Further, thresholds related 
to the antigen are also set to kq. On the other hand, in general, the antibody fa 
cannot recognise the antigen. Then, by assuming that the antigens enter the system 
at a rate Kj per unit time, the differential equation for the antigen A is given by 



dA 
~dt 



-K 1( r A (t)A + K 7 , 



(0-7) 



a A (t) = m 12 (t)fa(t) + m 13 (t)fa(t). 




(a)Case1 



(b)Case2 



Fig. 2. invasion of antigen 



Case 2 

An antigen A interacts only with 
the antibody fa. Then, the antibodies 
fa recognise the antigen, but antibod- 
ies fa and fa cannot recognise the anti- 
gen, in general. Let us define Kij as the 

threshold above which fa recognises fa. Here i and j take integer values between 1 
and 4 where i = 4 denotes the antigen. Then the equation for the antigen A is 



dA 

— = -K l( jA{t)A + K 7 , 

<y A {t) = 0(fa - K 4A )fa. 

The sensitivity 0, of the ith clone is modified as 



(0-8) 



(0-9) 



where k is the strength of the interaction between fa and A and is set to k = SA<5i,i 

50 for any i,j = 1,2,3 and ^4 = 



K 



3,1 



K 



with sa = 0.3. We assume Ki 
= «i = 10 for i = 1, 2, 3. 

We have performed numerical calculations in these two cases and have found the 
following results about the response of the system. As is shown in Ref. 9), there is 
a differentiating state with two L clones and one S, which responds to the antigen 
much better than the other states. In both cases 1 and 2, if the input rate Kj of the 
antigen is neither very small nor very large, the system modifies itself by evolving 
to that differentiating state in which the relaxation time of the antigen is shortest 
among differentiating states. This phenomenon is regarded as a kind of short term 
memory of the invasion by the antigen. However, in both cases, if if 7 is sufficiently 
large, the concentration of antibody which can interact with the antigen converges 
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to a fixed point, while the antigen concentration grows exponentially with the result 
that the system is destroyed by the antigens. The dependence on the initial condi- 
tion of this behaviour is less sensitive than that on K7. Further, we have found that 
the network is destroyed by the antigen more easily in the case 2 than in the case 1. 

Model without Threshold 

Next, we study the response of the system without the threshold to invasion by 
the antigen in cases 1 and 2. We introduce a strength parameter s of the connectivity, 
and set the elements of the connectivity matrix as 

my = sm% = s(l - 6 itj ). (0-10) 

As s increases from 0, the system changes as a fixed point — ► limit cycle 1 — > chaos 
— ► limit cycle 2 — > a fixed point. If the state is either limit cycle 2 or chaotic, the 
whole network is activated. On the other hand, in the limit cycle 1 state, a part of 
the network is activated, in the sense that one clone takes negligibly small values 
and does not affect the other clones. By fixing the initial value of the antigen and 
changing K7, we study the responses of these three states. In both Case 1 and Case 
2, for a small value of K7, the state of the system does not change very much, but 
for large Kj, each clone converges to a fixed point and the network is destroyed. The 
magnitude of the value of Ky that causes the breakdown of the system is of order 1 in 
Case 1, and is 1CP 1 in Case 2. In Case 1, for intermediate values of K-j, regardless of 
initial states, the attractor changes to a new periodic state(limit cycle 3) (Fig.||). This 
state is considered to be a better state against the invasion by the antigen because 
one of two clones which can interact with the antigen is in the long-period state, L. 
In Case 2, when K-j takes intermedi- 
ate values, the attractor changes to the 
chaotic state in which the clone 1 os- 
cillates in relatively longer pulse states. 
These behaviours occur for all initial 
states except for the state of the limit Fig. 3. Limit Cycle 3 

cycle 1 in which the clone 1 takes small 

values. Thus, as well as in Case 1, in Case 2 the system tends to take better config- 
urations to deal with the antigen. 

Therefore, whether there is the threshold or not, it turns out that in both cases 
1 and 2, the system changes to a desirable direction for moderate magnitudes of 
perturbation. Further, it turns out that the network is destroyed by the antigen 
more easily in Case 2 than in Case 1. The reason seems to be that the effect of the 
antigen on the network is more direct in Case 2 than in Case 1. 

Loosely Connected 3-clone Units 

Here, to study the effect of thresholds on the response to the invasion by antigens, 
we investigate one-dimensionally connected basic units, where each basic unit is 
composed of a 3-clone closed system. 
A. Basic units with threshold 
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As an initial state of a basic unit, we take a limit cycle in which the con- 
figuration of clones is (L, L, S). When two basic units are connected, for each 
unit the network architecture is similar to Case 2 of the invasion by the anti- 
gen. See Fig.0(a). Thus, from the above result in Case 2, if the interaction be- 
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Unit 1 Unit 2 Unit 3 



Unit 1 Unit 2 Unit 3 



(b) 

Fig. 4. Connection of 3-clone units 



(C) 




(1)ku=0 



(2)ku=50 



Fig. 5. K-j dependence of the magnitude of 
fluctuation in clone 3 and 8 relative to that 
in clone 



tween units is large, the combined system is destroyed when the amplitude in ei- 
ther unit becomes large. Hence, hereafter, we consider loosely connected systems. 
In the case of two units, the clone which 
is connected to the other unit tends to 
stay in the long pulse state, and the 
basin of attraction of this type of solu- 
tion becomes large. It turns out that 
this is a characteristic feature for the 
loosely connected systems. Now, let us 
study an open chain of three basic units: 
unit 1, unit 2 and unit 3. See Fig.||](b). 
Let su and k\j be the strength of the 
connection between units and that of 

the threshold between units, respectively. To see the response behaviour quanti- 
tatively, we have calculated in limit cycle states the differences between the phases 
of oscillations in the presence and the absence of the antigens. We have also changed 
Ato (the threshold between clones) and kjj- As a result, we have observed that the 
disturbance is reduced more in unit 3 than in unit 2. However, significant difference 
between in the cases of k\j = and hjj > has not been observed. See Fig. 5. Fur- 
ther, we have found that the disturbance is not always reduced as kq is increased, 
and that the threshold between clones provokes complicated dynamical behaviours. 
Therefore, to see the effect of the threshold more clearly, we investigate the system 
composed of the basic units without thresholds. 
B. Basic units without threshold 

We set Ko = 0. In this case, for s = 2.5 and sjj = 0.05, there exist limit cycle 
states in both cases of k\j = and 50. As in case A, setting Aq = 80, k\ = 10 and 
sa = 0.3, we have investigated the response behaviour of these limit cycles quan- 
titatively. We have observed that the disturbance again hardly spreads to unit 3. 
Further, we have found that the disturbance in unit 2 is reduced considerably for 
kjj = 50, while it is still large for kjj = 0. See Fig. 6. 



These results show that the thresholds have the effects of reducing the magni- 
tudes of interactions substantially, and of enhancing the independence of each unit. 
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(l)ku=o 



(2)ku=50 



Fig. 6. K-r dependence of the magnitude of 
fluctuation in clone 3 and 8 relative to that 
in clone 



Discussion 

In this letter, we have analysed the 
model constructed from the point of 
view of N. K. Jerne that antibodies look 
upon other antibodies as internal images 
of antigens. It has been pointed out as 
a big fault in his network theory that 
the stability of the network is not guar- 
anteed. There have been several stud- 
ies on this subjecti'iS. However, it 
seems that this problem has not yet been fully investigated. The result on immune 
responses obtained in the first part of this paper implies that the system does not 
always change randomly but can change in an ordered way in some cases. This may 
be a refutation of the above criticism of the network theory. Since this result has 
been obtained in systems with few degrees of freedom, it is necessary to study larger 
systems. This is now under investigation. 

Now, let us discuss the effects of thresholds of concentrations. The result in 
the second part of this paper shows that the structural stability of the network is 
increased by the existence of thresholds. It has been pointed out that a cross-linking 
structure made by antibody molecules is very important in order that phagocytes 
can catch antigens in real immune systems 0. It is natural to assume that the for- 
mation of the three-dimensional cross-linking structure by antibodies depends on the 
concentrations of antibodies and antigens. Thus, the introduction of the threshold 
of the concentrations not only gives a desirable effect which makes the network more 
stable, but also can be interpreted as taking one of the mechanisms of real immune 
responses into account theoretically. It would be very interesting to know whether 
the threshold provokes the localisation of the immune response in larger systems. 
This is a future problem. 



We are grateful to D. Broomhead and J. Huke for critical reading of the manuscript. 
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